Dimensional regularization applied to nuclear systems 
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We apply for the first time the dimensional regularization procedure to treat an ultraviolet diver- 
gence occurring in the framework of the nuclear many-body problem. We consider the second-order 
correction (beyond the mean-field approximation) to the equation of state of symmetric nuclear 
matter with a zero-range effective interaction. The unphysical ultraviolet divergence generated at 
second order by the zero range of the interaction is properly removed by the regularization technique 
and the regularized equation of state (mean-field -I- second-order contributions) is adjusted to a ref- 
erence equation of state. One of the main practical advantages of this procedure, with respect to a 
cutoff regularization, is to provide a unique set of parameters for the adjusted effective interaction. 
This occurs because the regularized second-order correction does not contain any cutoff dependence. 
The encouraging results found in this work indicate that such an elegant and rigorous technique to 
remove the divergent terms and to generate regularized effective interactions is likely to be applied 
in future in the framework of beyond mean-field models for finite nuclei. 
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Perturbation theories are employed in many domains 
of physics and regularization and renormalization tech- 
niques are adopted in cases where the inclusion of higher- 
order terms with respect to the leading contribution in 
the perturbative expansion generates divergences. This 
type of divergences are well known for instance in par- 
ticle physics in the context of quantum field theories 
[ij. In several frameworks of many-body physics, such 
as in nuclear physics and in atomic physics, other types 
of divergences occur in mean-field-based models like the 
Bogoliubov-de Gennes or the Hartree-Fock-Bogoliubov 
theories if a zero-range interaction is employed in the 
pairing channel to treat a superfluid many-fcrmion sys- 
tem Apart from this specific case related to the 
pairing interaction, in the perturbative treatment of the 
many-body problem an ultraviolet divergence always ap- 
pears if a zero-range interparticle interaction is used when 
higher-order terms are included beyond the mean-field 
approximation (which corresponds to the leading order in 
the perturbative solution of the Dyson equation). This 
unphysical divergence has been analyzed in a previous 
work (Hj. 

We develop here a technique which is currently adopted 
in the context of quantum field theories, the so-called di- 
mensional regularization, for the second-order correction 
beyond the mean-field approximation within a specific 
case of the nuclear many-body problem. As done in Ref. 

we consider the equation of state (EoS) of symmet- 
ric nuclear matter with a contact interaction gS{fi — 7^2) 
where the coupling constant g is density dependent and 
contains three parameters, to, and a: g{p) = to + ^ p"'. 
This corresponds to the — model of the nuclear ef- 
fective Skyrme interaction In Ref. Q we have al- 
ready analyzed in this specific model the nature of the 
divergence related to the second-order correction of the 
EoS. The divergent term has been calculated analytically 
and, by deriving its asymptotic expansion, it is possi- 



ble to show that it goes linearly with the momentum A 
which is introduced as a cutoff regulator. A cutoff regu- 
larization has been applied to absorb the divergence by 
means of a fit of parameters in the corrected mean- field -I- 
second-order EoS. In this work we treat the same diver- 
gence with a dimensional regularization technique which 
is for the first time applied to the nuclear many-body 
problem. The dimensional regularization technique is an 
elegant procedure which has the advantage of preserv- 
ing symmetry laws and of including high-energy effects 
which are sharply discarded with a cutoff regularization. 
This technique has been introduced in the framework of 
the electroweak theory 0-0] and consists in replacing the 
dimension of the divergent integrals with a continuous 
variable d. The main idea is that, if an integral diverges 
in a given integer dimension, the result may be finite 
by replacing the integer dimension with a non integer 
d. One then performs a kind of analytic continuation in 
the dimension to return to the initial integer value [lol |. 
The dimensional regularization eliminates power-law di- 
vergences, isolates logarithmic divergences and regular- 
izes infrared divergences [ll|. Together with the contin- 
uous variable d, a regulator e (which is dimensionless) is 
introduced so that, when e ^ 0, the dimension of the 
integral comes back to the initial integer value. In ad- 
dition, an auxiliary scale ^ is included to maintain the 
correct dimensions of the physical quantities which are 
calculated. After regularization, a renormalization is ap- 
plied by a minimal substraction procedure to remove the 
regulator e which appears as a 1/e pole if the divergence 
is logarithmic. All the physical observables are indepen- 
dent of the auxiliary scale fi due to the renormalization 
group equation fidS/dfj, = 0, where 5 is a generic observ- 
able. 

Let us consider the second-order correction beyond the 
mean- field EoS in our io~^3 model for symmetric nuclear 
matter Q . We write its generalized expression where we 
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introduce a continuous dimension d in the integral and the auxihary scale fi. In a box of volume n one has, 
I 
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In our case the regulator e can be written in terms of (i as e = 3 — d. When e — > 0, d returns to the integer value 3. 

I 

By making some manipulations and by using the In dimension d, all massless integrals are regularized 
Schwinger's proper time representation of Feynman in- to zero For instance, 

tegrals the following compact expression may be derived, 

A 



Cd{^i.p)g'{p) \ d'q Mq), (2) 
c, 



1= / d'^q d'^k2 ^ = 0. 
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where the domain of integration is 
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\ki\ < kp, \k2\ < fcF, \ki + g1, |fc2 - 9l > kp (3) 



By rewriting / as the sum of two integrals Ii and I2 
(by splitting the integration in two regions). 



with kp ~ (^^p) ■ The integrand in Eq. (2) can be 



written as 
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(4) the second-order correction reads 
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The quantity B is finite for \q\ < 2 when d ^ 3 and one A = , , ^''^,,^ T((j). By using hypergeometric func- 
can show that its value is equal to 
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The quantity A may be written as a function of d, 
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tions jl5l |. after straightforward but lengthy manipula- 



tions one can write 



T{d) 
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lim T(.). 16.^-1 + I log2). (9) 



T{d) converges for < d < 4 (with a pole at d = 4) as shown in Fig. 1. For d = 3 the divergence has 
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been removed by the regularization (the value of T(3) is 
positive and finite). 
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FIG. 1. T{d) as a function of the space-dimension d. 

By taking d = 3, the mean-field + second-order EoS 
finally reads 



E 



3^2 /37r 
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(10) 



One can notice that in Eq. (|T0|) the auxiliary scale /i 
(see Eq. (1)) has disappeared in the second-order correc- 
tion. This occurs because in our case the divergence is 
not logarithmic (no 1/e-poles are found after regulariza- 
tion). That means that we do not need to renormalize 
after regularization. We observe also that, as expected, 
the power-law divergent terms have been removed and 
do not appear in the final expression. Finally, if we 
compare Eq. (10) with the asymptotic expression of the 
second-order correction previously obtained with a cut- 
off regulator (that is, the asymptotic expression of the 
result shown in Ref. [11), we observe that the effect 
of applying the dimensional regularization technique is 
not only to remove the divergent term (the linear term 
in A is now absent), but also to renormalize the finite 
part: with the cutoff regularization, the numerical fac- 
tor multiplying 64'!r^C3{p)g'^{p) in the finite part is equal 
to jig ( — 11 + 21og2) whereas in Eq. (10) it is equal to 
jig (1 + 2 log 2). This difference comes out from the sub- 
straction procedure of the integrals Ii and I2 which renor- 
malizes also the finite part. Notice that this substraction 
must be done carefully and the domain of integration 
must be taken correctly as the intersection of C/ and 
|g| > 2 (|g| < 2) for h {h)- 

In the upper panel of Fig. 2 one can see the regu- 
larized second-order EoS (dashed line) with respect to a 
reference Skyrme mean- field EoS (full line). As in Ref. 
0, the SkP 14 1 mean-field EoS has been chosen as a 
benchmark. The regularized second-order EoS has been 
calculated with the SkP parameters. We notice that the 
second-order correction at the saturation point is ^ - 8 
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MeV, which is comparable to what obtained in Ref. in 
the case of a cutoff calculation for a value of the momen- 
tum cutoff A = 1.5 fm^^. We also notice that the profiles 
of the two second-order corrections (with a cutoff calcu- 
lation and with the dimensional regularization) are quite 
different. The removal of the divergent part in the latter 
case leads to an equilibrium point located at 0.167 fm~'^ 
which is close to the saturation point of matter. This was 
not the case for the cutoff calculation where the minimum 
was appreciably shifted by the second-order term. The 
shift was actually due to the density dependence of the 
integral which multiplied the factors C3 and g^. In the 
present case the density dependence of the second-order 
correction is contained only in the coefficients C3 and g^ 
and does not modify in practice the position of the equi- 
librium point which is provided for symmetric matter by 
the first-order contribution. 

Since SkP is a phenomcnological interaction where the 
parameters have been adjusted at the mean-field level, we 
have to readjust the parameters for the second-order and 
regularized EoS to have a reasonable curve for the en- 
ergy. We choose as a reference for this fit the mean-field 
SkP-EoS as done in Ref. [![. Upon fitting the regular- 
ized EoS, a unique set of Skyrme parameters (no cutoff 
dependence) is generated (parameters are listed in Ta- 
ble 1). The fit is performed on 15 points on the SkP 
mean- field curve up to a density of 0.30 fm"'^. The pro- 
cedure is a minimization with the following definition 
of X^, = E^iiEmedii) - Eref{i))/AE,]\ The errors 
AEi in the denominator have been chosen equal to 1% 
of the reference SkP energies. The resulting per point 
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TABLE I. Refitted parameters compared with the SkP values 
for to, t-i and a. 



to (MeV fm^') 


t3 (MeV fm3+3") 


a 


Skp -2931.70 


18708.97 


0.16667 


New parameters -2893.730 


17977.603 


0.131 



is equal to 0.547 indicating that, on average, the fitted 
points deviate less than 1% from the reference SkP-EoS. 
We can conclude that the fit is of very good quality. This 
value of the P^r point cannot be compared with the 
corresponding numbers found in Ref. Q in the case of 
the cutoff regularization because a different expression of 



has been used in the two works. 



FIG. 2. Upper panel: Regularized second-order (dashed line) 
and mean-field SkP (full line) EoS. Lower panel: Regularized 
and refitted second-order EoS (dashed line) compared with 
the SkP mean-field EoS (full line). The star indicates the 
position of the saturation point. 
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The second-order regularized and refitted EoS is shown 
in the lower panel of Fig. 2. The saturation point is 
located at 0.161 fm""^. The curve is practically super- 
posed to the mean- field SkP-EoS. The resulting pres- 
sure P(p) ~ P^-j^ [^(p)] ^^'^ incompressibility modulus 

Koo{p) = 9p^^ ^1^° shown in Figs. 3 and 4, 

respectively. We stress that these quantities do not enter 
explicitly in the fit. The agreement with the correspond- 
ing SkP values is good except at densities larger than ^ 
0.25 fm~'^. The region around the saturation point has 
thus reasonable properties. Based on these encouraging 
results, work is presently in progress to extend the ap- 
plication of the dimensional regularization procedure to 
the whole Skyrme effective interaction by including also 
the velocity-dependent terms which are more strongly 
divergent [16|. This opens new prospectives for future 



applications of dimensional regularized Skyrmc-type ef- 
fective interactions adjusted at a beyond mean-field level. 
These interactions would be well adapted to be used in 
the framework of beyond mean-field approaches for finite 
nuclei. 



FIG. 3. Pressure calculated at second order with the refitted 
parameters (dashed line) compared with the SkP-mean-field 
curve (full line). 
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FIG. 4. Incompressibility modulus calculated at second or- 
der with the refitted parameters (dashed line) compared with 
SkP-mean-field curve (full line). 
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